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ABSTRACT 

A box model is used to study first order Fermi acceleration driven by magnetic re- 
connection. It is shown, at least in this simple model, that the spectral index of the 
accelerated particles is related to the total compression in the same way as in diffusive 
shock acceleration and is not, as has been suggested, a universal E~ 5 / 2 spectrum. The 
acceleration time-scale is estimated and some comments made about the applicability 
of the process. 
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1 INTRODUCTION 

It has been suggested that the convergent flows as- 
sociated with magnetic reconnection may drive a pro- 
cess of first order Fermi accel e ration . In particular, in 
" Ide Gouveia dal Pino fc LazarianI (|2005f ) it is argued that 
such a process can occur and naturally produce an E~ 5 ^ 2 
power-law differen tial energy spect rum of accelerated parti- 
cles. More recently I Gianniosl (2010) has invoked this process 
, to explain the acceleration of the UHECRs (co smic rays with 
■ E > 1 19 e V). It has also b een c onsidered in ICerutti et al.l 
(2011) and ILazarian et alJ (|201 lh . It is thus of interest to 
examine this idea in more detail. 



2 THE MODEL OF DE GOUVEIA DAL PINO 
AND LAZARIAN 

The basic idea is quite simple. The reconnection layer is ap- 
proximated as a thin planar sheet with plasma flowing in 
from both sides at the reconnection velocity Kec- Charged 
particles with gyro-radii much larger than the thickness of 
the sheet are then assumed to freely diffuse in the inflow- 
ing plasma and, much as in the standard theory of diffusive 
shock acceleration, to gain energy every time they cross the 
reconnection layer. As in the theory of shock acceleration, 
this relies on using a mixed system of coordinates where par- 
ticle momenta (and associated energies) are measured with 
respect to the local co-moving plasma frame. This has the 
advantage that all magneto-static scattering effects in the 
plasma on either side of the layer, responsible for isotropis- 
ing the particle distributions and producing spatial diffu- 



sion, are then energy-neutral and the energy changes are 
entirely due to the frame-change on crossing. It should be 
noted that as long as the distributions are close to isotropy 
and the sheet is assumed to be thin, direct acceleration by 
the electric fi eld in the reconnection layer as suggested by 
ISpeiserl (|l965t) is unimportant except for a very small set of 
particles on Speiser orbits. 

An argument essentially identical to that given by Bell 

(1 19781 ) for the case of shock acceleration then shows that on 
each crossing of the reconnection sheet the particle momen- 
tum or energy increases, averaging over an isotropic distri- 
bution, by an amount 

Ap = A£ = 4K 2£ 

p E 3 c U 

or twice this for a cycle of two crossings (bringing the particle 
back to the side it started on). The fractional gain in momen- 
tum (Ap/p) resulting from the frame change on crossing the 
discontinuity is easily found to be 1 + fi/3 where f3 = 2V rec /c 
is the dimensionless velocity jump across the discontinuity 
and fi = cos 6 is the cosine of the angle at which the par- 
ticle crosses with respect to the normal. Averaging over an 
isotropic distribution, (fj,) = 2 ji 2 d^i — 2/3 then gives the 
above result. 

de Gouveia Dal Pino and Lazarian then argue that "the 
escape probability of a particle from the acceleration zone is 
similar to the one computed for a shock front" and that the 
probability (per cycle of crossings) of staying within recon- 
nection region is 
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Figure 1. A cartoon sketch of a generic magnetic reconnection 
system 



Now if the energy per cycle increases by a factor 1 + a and 
the probability of rem a ining in the system per cycle is 1 — ft 
it is easy to show lBeill (1 19781 ) that the steady state differen- 
tial energy spectrum is N(E) oc £ ,_1- ' 3 / a from which they 
deduce a spectrum of E~ s ^ 2 . 

There are two serious objections to this analysis. First, 
in calculating the energy gain they assume that the flow is 
entirely compressive and neglect the outflow from the recon- 
nection region which must be associated with a divergence of 
the flow field (and associated energy losses). Secondly, their 
assumption that the escape rate is the same as that from a 
shock is entirely without justification. Even it it is accepted 
as an order of magnitude estimate, it certainly does not jus- 
tify the precise value required to get the spectral index of 
—5/ 2. Nevertheless the idea is interesting and worth pursu- 
ing. iGianniosl (|2010) considers the ultra-relativistic version 
of the process, but focuses on the maximum attainable en- 
ergy and does not consider the spectrum. 



3 A BOX MODEL TREATMENT 

Let us model the reconnection process as a 'black box' sys- 
tem into which highly magnetized plasma flows symmet- 
rically from above and below and out of the sides of which 
emerge two 'jets' of plasma with little internal magnetic field 
(see Figure [l}. 

There are two conservation laws that can be applied to 
constrain this system. First, the flux of mass in has to equal 
the flux out if the system is quasi-steady, 



piUiAi = p 2 U 2 A 2 = M, 



(3) 



where the incoming plasma has density pi, velocity Ui and 
the cross-sectional area of the reconnection region for the 
inflow is Ai and similarly for the outflow. 

Secondly, the energy flux into the system has to equal 
that out. Noting that the magnetic energy density is related 
to the local Alfven speed, Va by 



B 1 T/2 
2i^ = 2 pV *> 



(4) 



and if we assume that magnetic energy density dominates 
in the inflow and kinetic energy density in the outflow, this 
gives 
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(5) 



from which the well know result follows that the outflow 
velocity is of order the Alfven speed in the inflow. 

Let us now consider the acceleration of particles due to 
the convergence of the flow in the reconnection region (which 
we denote as R) following the method of the si mplified 'box 
model ' of shock acceleration as discussed in iDrurv et al.l 
The upwards flux of particles in momentum is simply 
given by 



$ = 



W/(p)§("V -V) 



(6) 



and if we assume that for sufficiently high-energy particles 
the particle distribution is essentially uniform within the 
reconnection region then we get 



$ = 



4np s 



f(p) / (-V • U) 
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f( P ) [2A x Ui - 2A 2 U 2 ] 



(7) 
(8) 
(9) 



As in the shock case the steady spectrum can now be found 
by equating the divergence of the upwards flux in momen- 
tum to the losses carried out by the outflow, 



^ = -2A 2 U 2 4irp 2 f{p) 

which, after a little algebra, gives the result 
9 In/ _ 3AiUi 



9 hip 
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(10) 



(11) 



In fact this is precisely the same result as in shock accel- 
eration if we express it in terms of the compression in the 
system. Using the fact that 



A1U1 = M/pi 



A2U2 = M/p 2 



(12) 



(13) 



one easily finds that 

9 In/ 3r _ p2 

dlnp r — 1' pi 

The question is thus how compressive can the process of 
magnetic reconnection be? If little energy is used to heat the 
plasma and the conversion is essentially one of magnetic en- 
ergy to kinetic energy it is clear that it can be very compres- 
sive. The only constraint appears to be that the outflow not 
be significantly over pressured relative to the environment, 
and for a strongly magnetized inflow this is a very weak con- 
straint. Thus in general one expects rather large values of 
the compression ratio, certainly larger than the four canon- 
ically assumed for shocks, and thus extremely hard spectra 
tending in the limit to 



f(p) oc p 



(14) 



or a differential energy spectrum of N(E) oc E^ 1 . Of course 
such hard spectra rapidly lead of energy divergences and 
non-linear effects, but at least at the test-particle level it 
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would appear that reconnection regions may be efficient sites 
for producing hard spectra by Fermi acceleration. 

It is also quite straightforward to derive the acceleration 
time-scale. If the accelerated particles penetrate a distance 
L\ 'upstream' into the inflow and L2 'downstream' into the 
outflow then it is easy to show that the acceleration time 
scale is 

_ 1} A l L 1 + A2L2 + W 
* acc _ 6 A 1 U 1 -A 2 U2 1 ' 

where W is the volume of the reconnection region; this again 
is a straightforward generalization of the result for shock 
acceleration and physically just expresses the fact that the 
upward flux in momentum takes a certain amount of time 
to fill the acceleration volume. 

The above analysis, as is the case also in most discus- 
sion of the test-particle theory of shock acceleration, glosses 
over the question of where the particles come from (the in- 
jection problem). It simply assumes that there is some in- 
jection process operating at energies lower than those being 
discussed which produces a population of supra-thermal par- 
ticles which can then be further accelerated by this Fermi 
process. In the case of reconnection this could be plasma 
physical processes related to the micro-physics of the ac- 
tual reconnection, or it could be a pre-existing population 
of ambient particles (ambient cosmic rays). 



theory possible by allowing the pairs to annihilate in the 
compressed reconnection sheet. This may be relevant to ac- 
celeration by reconnection in pulsar winds. In this case an 
extension of the model to include synchrotron losses would 
be required. 
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4 CONCLUSIONS 

This approximate analysis suggests that reconnection can 
indeed drive a process of Fermi accel eration in clos e 
analogy to diffusive shock acceleration dKrvmskii Il97' 



iBlandford fc Ostriker! Il97l I Axford et al. 1 119771 : : iBell 1 19781 ) 
even to the extent of having the same relation between total 
compression and the spectral index of the accelerated par- 
ticles (at least to a first approximation). Because there is 
no constraint from momentum conservation (the total mo- 
mentum being identically zero by symmetry) the total com- 
pression will normally be quite high leading to rather hard 
spectr a. The suggestion bv lde Gouveia dal Pino fc Lazarian! 
(2005) of a universal E~ 5 ^ 2 spectrum is not supported by 
this analysis. 

The acceleration time scale is similar to that for shock 
acceleration, and the constraints on maximum energy due 
to the finite age and size of the reconnection region will 
thus be comparable to those in shock acceleration. In gen- 
eral the much smaller scale of reconnection regions compared 
to shock structures will thus lead to significantly lower max- 
imum energies (the dominant limiting loss process being dif- 
fusion out of the sides of the reconnection region). Strictly 
the model developed here requires that the particle diffusion 
length scale in the inflow be large compared to the thickness 
of the reconnection sheet, but small compared to its lateral 
extent; thus it may only apply in a rather restricted energy 
range (this should be contrasted with the case of shock ac- 
celeration where the scale separation is much more extreme) . 
However the particles accelerated in the reconnection regions 
may be available as a seed population for further shock ac- 
celeration in many systems, and the fast time-scales and 
hard spectra may be of interest in models of flaring activity. 

Finally it is interesting to note that for reconnection in 
an electron-positron pair plasma infinite compression is in 



